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Abstract 

We analyze the big bounce transition of the quantum FRW model in the setting of the non- 
standard loop quantum cosmology (LQC). Elementary observables are used to quantize composite 
observables. The spectrum of the energy density operator is bounded and continuous. The spec- 
trum of the volume operator is bounded from below and discrete. It has equally distant levels 
defining a quantum of the volume. The discreteness may imply a foamy structure of spacetime 
at semiclassical level which may be detected in astro-cosmo observations. The nonstandard LQC 
method has a free parameter that should be fixed in some way to specify the big bounce transition. 
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I. INTRODUCTION 



It is commonly expected that the cosmological singularity (CS) problem may be 
resolved in a theory which unifies gravity and quantum physics. It seems that recent de- 
velopments concerning quantization of cosmological models by making use of loop geometry 
may bring solution to the problem. It consists in turning the classical hig bang singularity 
into hig bounce, BB, transition. There exist two methods to address the issue: standard loop 
quantum cosmology (LQC) and nonstandard LQC The former has been developed during 
the last decade (see p, m ^'^^ references therein) and has been inspired by the loop quan- 
tum gravity, LQG (see [7H9| and references therein). The latter has been proposed recently 
10-13| and seems to be related to the reduced phase space quantization of LQG fl3 |. 

The standard LQC means basically the Dirac method of quantization, which begins with 
quantization of the kinematical phase space followed by imposition of constraints of the 
gravitational system in the form of operators at the quantum level. Finding kernels of these 
operators helps to define the physical Hilbert space. In the nonstandard LQC (our method) 
one first solves all the constraints at the classical level to identify the physical phase space. 
Next, one identifies the algebra of elementary observables (in the physical phase space) and 
finds its representation. Then, composite observables are expressed in terms of elementary 
ones and quantized. Final goal is finding spectra of composite observables which are used to 
examine the nature of the BB phase in the evolution of the universe. 

In what follows we restrict our considerations to the fiat Friedmann- Robert son- Walker 
(FRW) model with massless scalar field. This model of the universe includes the initial 
cosmological singularity and has been intensively studied recently within the standard LQC. 

In our recent paper [lo[ we have argued that the resolution of the singularity offered by 
LQC requires specific value of the fundamental length. The size of this length has not been 
determined satisfactory yet. Present paper is an extended version of [lH and together with 
the classical formalism [l2| specifies our nonstandard LQC (when applied to FRW model). 

In order to have our paper self-contained, we recall in Sec. II some aspects of the classical 
formalism of our nonstandard LQC [l2[. In Sec. Ill we present our quantization procedure. 
It consists in finding representation for composite observables like the energy density and 
the volume function, and calculating their spectra. We conclude in the last section. 

In the Appendix, we give an extended motivation for our paper. In particular, we justify 
the need for quantization of our loop FRW model [12], which is free from a cosmological 
singularity already at the classical level due to the modification of the classical Hamiltonian. 

II. CLASSICAL LEVEL 
A. Hamiltonian 

The gravitational part of the classical Hamiltonian, Hg, of the fiat FRW model may be 
presented in the form 121 

i/,= limiff, (1) 



where 



^^'^ = -£0jr^T.^''''Tr[h\^^h^\h^Y\h^Y'hi-\ih^^^^^^ (2) 

ijk 
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and where V = \p\2 = a^Vo is the volume of the elementary cell V C S (E is spacelike hyper- 
surface); j9 is a canonical variable; the metric of = FRW model is: ds"^ = — iV^(t) dt"^ + 
a^(t) {dx"^ + dy"^ + dz"^), where a is the scale factor and N denotes the lapse function; e*-^^ 
is the alternating tensor; 7 is the Barbero-Immirzi parameter, and the holonomy function 
h^il^\ along straight line of coordinate length proportional to ;u/|a|, reads 

hi^^ = cos(/i/3/2) I + 2 sin(/i/3/2) r^, (3) 

where Tk = —iak/2 {ak are the Pauli spin matrices). 

The total Hamiltonian for FRW universe with a massless scalar field, 0, is found to be 

H = H^ + H^, (4) 



where = A^p||p|~5/2, and where and are canonical variables. Equation (j4]) satisfies 
the condition if ^ as it corresponds to the scalar constraint of general relativity; the 
Gauss and diffeomorphism constraints are equal zero (in the strong sense) due to the choice 
of gauges. 

Making use of ([3]) we calculate ([2]) and get the modified total Hamiltonian Hg'^^ corre- 
sponding to (jl]) in the form T3 

hW/N- 3 sin^(A/3) ^ 

where 

/3:=^^ -=br (6) 

are the canonical variables of the so called improved scheme, and where c and p are canonical 
variables. The variable /5 = 70/0 corresponds to the Hubble parameter d/a, and v^/^ = 
aV^^^ is proportional to the scale factor a. The relationship between the coordinate length 
/i (which depends on p) and the physical length A (which is a constant) reads: X = fx = 
/X a \/q . 

At this stage, it should be emphasized that (|5]) presents a purely classical Hamiltonian 
modified by the holonomy (|3]), i.e. Eq. ([5]) includes no quantum physics. Contrary, in the 
standard LQC ([5]) is called a semi-classical or effective Hamiltonian and is interpreted to 
include some 'imprints' of quantization. In our nonstandard LQC quantum physics enters 
the formalism only when quantizing the algebra of observables (presented in the subsequent 
section). 

The complete Poisson bracket for the canonical variables (/3,f,0,P(^) is defined to be 



{■,-}:=47rG7 



d- d- d- d- 



d(f) dp^ dp^ 



d- d- d- d- 

d(3 dv dv (9/3 

The dynamics of a canonical variable ^ is defined by 

e:={e,/7W}, ^e{f3,v,(P,p4, (8) 

where ^ := d^/dr, and where r is an evolution parameter. 

The dynamics in the physical phase space, J^p'^y^, is defined by solutions to (IHl) satisfying 
the constraint H^^"* 0. The solutions of ([8]) ignoring the constraint are in the kinematical 
phase space. 
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B. Elementary observables 

A function, O, defined on phase space is a Dirac observable if 

{O,H^^^}^0. (9) 



Equation ([9]) leads to [12 



^— — - . cos(A/3) — ^ — = 0, (10) 

where := AnG/S. Solution to (|TOi) in J-'^^^^ is found to be H 



Oi:=P^, O2:=0-^^^arth(cos(A/3)). (11) 
One may parameterize J-'^^^ by the elementary observables (9i and O2, and define the 



Poisson structure as follows [12 



For simplicity we use the same notation for the Poisson bracket in (11211 and ([7]). 



C. Compound observables 

In what follows we consider functions on physical phase space like the energy density 
and the volume, which describe singularity aspects of our cosmological model. Considered 
functions are expressed in terms of elementary observables and an evolution parameter 0. 
They become observables for each fixed value of 0, since in such case they are only functions 
of observables [l2| . 

An expression for the energy density p of the scalar field (f) reads jl2| 

p((f),X) = - 1 . (13) 

' 2 {KjXy cosh" 3k{(P - O2) ^ ' 

The volume, a geometrical function, is found to be 13 

v{<p, A) = K7A \Oi\ cosh3K(0 - C2). (14) 

III. QUANTIZATION 

By quantization we mean: (i) finding a self-adjoint representation of composite observ- 
ables, and (ii) calculating spectra of operators corresponding to the composite observables. 
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A. Representation of elementary observables 

In what follows we use two representations of the classical algebra f|T2l) . Namely, 

Oi ai/(x) := -I hdJix), O2 ^ d2f{x) := xf{x) := xf{x), (15) 

which leads to [Oi, 02] = —ihl, and 

O, difix) := xf{x) := xf{x), O2 02/(x) := -ihdj{x), (16) 

that gives [Oi,02\ = ihl, where x G M. Due to the Stone— von Neumann theorem all 
self-adjoint representations of the algebra f|T2l) are unitarily equivalent to the representation 
f fTSjl or f lT6|l defined on a suitable dense subspace of L^(R). In this sense our choice of 
representation for (|T2l) is unique. 



B. Energy density 

The representation (fT6|) is essentially self-adjoint on on the dense subspace D of the 
Hilbert space L^[— r, r], where r G M+, defined to be 

D:={fe C°°[-r,r] | f-\-r) = /(")(r), n G {0} U N}, (17) 
where /"^"^ := d^f/dx"'. The eigenvalue problem, C2/P = P fp, has the solution 

fp{x) = {2r)-^/^exp{ixp/h), p{k) := 2nhk/r, k e Z. (18) 

The spacing □ of neighboring eigenvalues 

D:=p{k + l)-p{k) = 27ih/r (19) 

can be made as small as desired by making r sufficiently large. Thus, one may say that the 
spectrum of O2 is continuous. 

In the representation (fT6|) the energy density operator reads 

p — ^ ^ (20) 

2 (k7A)^ cosh^ 3k(0 + i hd^) 

Since C2 is essentially self-adjoint on J> := {fp(k)}k&i we may apply the spectral theorem 
to get 

pfp = p{<P,^,p)fp, (21) 
where ^ ^ 

'^'^ 2 (k7A)2 cosh^ 3k(0 — p) ' ^ 

and where p(0, A,j9) is the eigenvalue corresponding to the eigenvector fp. 

Our results mean that classical ( IT3i) and quantum (l22l) expressions for the energy density 
coincide. It is clear that the maximum density, pmax, reads 

Starting from the representation fflSll . instead of f|T6l) . we would get the quantum model of 
the energy density presented in [ll|, which is equivalent to the present one. 



5 



C. Volume operator 

To define the quantum operator corresponding to f , we use 

v = \w\, w:= K-fXOi cosh 3K{(f)- O2). (24) 

Thus, quantization of v reduces to the quantization problem of w. Quantization of the latter 
may be done in a standard way as follows 

wf{x) := K7A^ (^di cosh 3/t(0- 62) + cosh 3k(0- O2) Oi^fix), (25) 

where / G L^(R). For the elementary observables Oi and O2 we use the representation (|T5ll . 
An explicit form of the operator w is 

w f{x) = i — ^ — ^2 cosh3ft(0 ~ ^) S/tsinh 3fi;(0 — x)^ f{x). (26) 

Considerations simplify if we take / in the form 

f{x) := A e*^(^') cosh-^/2 3k(0 - x), (27) 
where ^4 G M, and where h is a. real-valued function. 

1. Eigenvalue problem 

Let us consider the eigenvalue problem for the operator w in the set of functions of the 
form ([27]). We have 

wf(x) = -K-^Xhcosh3K((f) - x)^^^^ f(x) =: bf(x), (28) 

ax 

where 6 G M is the eigenvalue of w. One may verify that a general form of h satisfying f l28p 
is given by 

2h 

h(x) = , arctan e^'^^"") . (29) 
Thus, a normalized satisfying (!28|) reads 

— exp (i^r^Kr arctan e^'*^'^"^) 
cosh*3K(^-.) ■ 



2. Orthogonality 
Making use of ( 130|) gives 



" vr cosh 3/t(0 — x) 
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We introduce a new variable y = e^'^^'^ and have 

9 /•(» exp I i ^o'fl arctan-u ) 
{/.!/.) = -/ " t . -^dy. (32) 

Another substitution tan z = y leads to 



TT 



It is clear that {fb\fa) = iff 

a — b = GK^'jXhm = SrrG'jXhm, m E (34) 

Thus, the set J-";, := { | a = 6 + STcG'-fXhrn; m G Z; 6 G M } is orthonormal. Each 
subspace C spans a pre-Hilbert space. The completion of each span Th-, V6 G M, 

in the norm of L^(M) gives L^(R). 

5. Self-adjointness 

Since {fb\wfa) — {wfb\fa) = {a — b){fb\fa), the operator w is symmetric on for any 
6 G M, because {fb\fa) = for a 7^ 6 due to the orthogonality of the set J^b- 

To examine the self-adjointness of the unbounded operator w, we first identify the defi- 



ciency subspaces, /C±, of this operator [15|, [16 



/C± := {^?± G I {g±\{w ± a)/,) = 0, V/„, G (35) 

where Db{w) := span J^, and where Db{w*) := {/ G L^(]R) : 3!/* = {f\wg), Wg G 

A(^)}. 

For each fa G Db{w) C L^(M) we have 



0= ((7±|(t&±^I)/,) = (a±z) / dxg±{x)Ux) =^ = Q = g_. (36) 

^—00 

Thus, the deficiency indices n± := dim[)C±] of w satisfy the relation: = = n_, which 
proves that the operator w is essentially self-adjoint on Db{w). 

4- Spectrum 



Due to the spectral theorem on self-adjoint operators [15|, ll6|, we may carry out quanti- 
zation of the volume as follows 

V = \w\ — y vfa := \a\fa. (37) 

A common feature of all J-f, is the existence of the minimum gap A := SnG'yhX in the 
spectrum, which defines a quantum of the volume. Let us discuss this issue in more detail. 
Let us denote the minimum eigenvalue of v by fmm- One can verify that Vmin = min{6. A— &}, 
where b G [0, A[. The spectrum consists of the union of {vmin + nA} and {—Vmin + {n + l)A}, 
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where n = 0, 1, . . . Only in two cases these two subsets are identical, namely when Vmin = 
or Vmin = , for which the minimum gap A is a constant gap between any two adjacent 
levels of the spectrum. Otherwise, the gap equals either A — 2vmin or 2vmin , and the 
minimum gap is the smaller one. One can verify that the case of any 6 G M reduces to the 
above case. 

In the limit A — ?• 0, corresponding to the classical FRW model without the loop geometry 
modification, there is no quantum of the volume. 

It results from flM|) that for 6 = and m = the minimum eigenvalue of v equals zero. 
It is a special case that corresponds to the classical situation when v = 0, which due to 
(|5]) means that = so there is no classical dynamics (for more details see [l2l)- Thus, 
we have a direct correspondence between classical and quantum levels corresponding to this 
very special state. It is clear that all other states describe bouncing dynamics. 

As the universe expands, a discrete spectrum of the volume operator may favour a foamy 
like structure which should finally turn into a continuous spacetime. The quantum of a 
volume may be used as a measure of a size, Xf, of a spacetime foam. One may speculate 

1 /3 

that Xf := _ [8TTG-fhX) ' . Thus, an astro-cosmo data that determine a size of 

spacetime granularity may fix the minimum length parameter A of LQC. That would enable 
making an estimate of the critical matter density pmax = l/2(fi:7A)^ corresponding to the 
Big Bounce. 

The granularity of volume should lead to the granularity of energy of physical fields. 
We suggest, making use of the de Broglie relation, that a specific particle representing a 
quantum of energy may have a momentum pi corresponding to its wavelength Aj such that 
Pi Xi = h. The detection of an ultrahigh energy particle with specific Pi may be used to 
determine Aj, and consequently set the upper limit for the fundamental length Xf. The set 
of parameters A, (for a set of particles) may be treated further as multiplicities of A/ in 
which case the greatest common divisor of all Xi would set the lowest upper limit for A/. 

5. Evolution 

The relation between eigenvectors corresponding to the same eigenvalue for different 
values of the parameter reads: 



f^+^ = ^^d,f^ ^ e-4o.f^_ (38) 



One may verify that 



w{(P + ^) = cosh i3K^)wi(P) + ^^^^M^9^^(0), (39) 



thus 



(AV(0 + ^)/a") = (AV(0)/.^)cosh(3«:^) + "-^^^^^{ft\dM<l>)ft) 
= a cosh {Sk^P) 6a, + {b- af^^^^^^^if^mt)- (40) 
Finally, an evolution of the expectation value of the operator w is found to be 

{f{(P)\w{(P + ^/^)/(0)) = A cosh 3k(V^ + B), (41) 
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where / := ^Oa/a, fa ^ J^b- One may verify that 



1 . X + Y 



A = sgn(X)VX2-y^ B = —\n — ^ , (42) 

OK X — Y 

where 

X:=)^K|a, r:=)^ (^m - l)(2m + 1) ' ^'^^ 

a a, m \ / \ / 

and where 6 = a + Gn^'yXh, 6 G M, m G Z, and |X| > 

We can see that the evolution of the expectation value of the operator w coincides with 
the classical expression (!24l) . 



D. Energy density operator in the basis J'h 

The operator p may be expressed in terms of the basis vectors from J^b. One may verify 
that 

(/a+mA|p/a) = < 4(;^ if TU = 0] (44) 

otherwise 
which leads to ^ ^ ^ 

" 4(«:7A)2^'^ ~ 8(/s:7A)2^'^-'^ ~ 8(/€7A)2^'^+'^' ^'^^^ 

where A := SnG'yhX . Therefore, p is bounded and does not commute with v. The latter 
is consistent with the Poisson bracket relation for these observables, which reads 

{m;,p} = 3fi;27Asinh3fi:(0- C2)p. (46) 

Therefore, the operators p and v cannot have common eigenfunctions so they provide two 
alternative ways of describing the system, either in terms of the eigenfunctions of p or v. 



IV. SUMMARY AND CONCLUSIONS 

Turning the big bang into the big bounce (BB) in the FRW universe is due to the 
modification of the model at the classical level by making use of the loop geometry. The 
modification is parameterized by a continuous parameter A. 

Each value of A specifies the critical energy density of the scalar field corresponding to the 
BB. The spectrum of the energy density operator is bounded and continuous. Classical and 
quantum expressions for the minimum of the energy density coincide. The spectrum of the 
volume operator, parameterized by A, is bounded from below and discrete. The expectation 
value of the volume operator coincide with the classical expression. The results concerning 
the volume operator may be extended to the area and the length operators [l2| . 

An evolution parameter, 0, does not belong to the physical phase space of our nonstandard 
LQC. Thus, it stays classical during the quantization process [lH. In the standard LQC, 
contrary to our method, is a phase space variable so it must be quantized [13] ■ Being a 
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quantum variable it may fluctuate so its use as an evolution parameter at the quantum level 
has poor interpretation. 

The Hamiltonin constraints are treated differently in the standard and nonstandard LQC 
at the classical level. In the former case the Gauss and diffeomorphism constraints are 
turned into zero by a suitable choice of gauges; the scalar constraint is solved only at the 
quantum level. In the latter case all the constraints are classically solved leading to the 
physical phase space. Thus, considered observables are not kinematical, but physical. This 
is why an evolution of geometrical functions like the volume may be used for testing the 
singularity aspects of a given cosmological model. The standard and our LQC methods 
give similar results in the sense that they lead to the Big Bounce transition. However, our 
method is fully controlled analytically as it does not require any numerical work. 

We believe that our nonstandard LQC may be related with the reduced phase space 
quantization of Loop Quantum Gravity (LQG) (see \l4^ and references therein). Finding 
the correspondence may help deriving LQC from LQG. 

There exist results concerning the spectrum of the volume operator obtained within LQG 
(see, e.g. [H, [l9[), but cannot be compared easily with our results due to the lack of a direct 
correspondence between LQG and LQC models. 



As there is no specific choice of A, the BB may occur at any low and high densities [10 
Finding specific value of the parameter A (and energy scale specific to the BB) is an open 
problem. It may happen, that the value of the parameter A cannot be determined, for some 
reason, theoretically. Fortunately, there is a rapidly growing number of data coming from 
observational cosmology that may be useful in this context. In particular, the detection of 
the primordial gravitational waves created at the BB phase may bring valuable information 
about this phase [20-22||. There exists speculation that the foamy structure of spacetime 



may lead to the dependence of the velocity of a photon on its energy. Such dependance is 
weak, but may sum up to give a measurable effect in the case of photons travelling over 
cosmological distances across the Universe [23|. Presently, available data suggest that such 
dispersion effects do not occur up to the energy scale 5 x 10^^ GeV [2^ so such effects may 
be present, but at higher energies. 

Various forms of discreteness of spacetime underly many approaches in fundamental 
physics. Just to name a few: noncommutative geometry 25l |. causal sets approach (26j . 



gravitational Wilson loops [27[, Regg e calculus |28[, path integral over geometries [29[, spin 



foam model [30[, and categories [31|. The discreteness may translate at the semi-classical 
level into a foamy structure of space. Such expected property of spacetime creates large 
activity in observational astrophysics and cosmology (see, e.g. Lorentz and CPT violation 



32 1 , dispersion of cosmic photons [3J] , electrons j33| and neutrinos j35| , birefringence effects 
36(1). Our results concerning the physics of geometry at short distances give some support 

to these approaches and expectations. 

For more concluding remarks, connected with motivation for our approach, we recommend 

reading the appendix section. 
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Appendix A: Motivation 

In what follows, we give an extended motivation for our paper. In particular, we justify 
the need for quantization of our loop FRW model [l2[, which is free from the cosmological 
singularity already at the classical level due to the modification of the classical Hamiltonian: 

If we intend to quantize canonically a system with Hamiltonian constraint, we can apply 
two methods: (i) Dirac's quantization - quantize kinematics ignoring the constraint, than 
impose the constraint at the quantum level (scheme D), or (ii) reduced phase space quanti- 
zation - solve the constraint at classical level, than quantize the system free of the constraint 
(scheme R). 

Suppose there is a theoretical reason to quantize a classical system with the Hamilto- 
nian constraint, but experimental/observational data with expected quantum effects are not 
available yet. This is the situation specific for gravitational systems. Theoretical reasons to 
quantize gravity are, e.g., singularities of classical theory and belief that all interactions can 
be unified. Which method should we choose to achieve the goal? Since the quantum data 
are not available yet, the best strategy seems to be applying both methods to compare the 
results. An agreement of the results would prove that the procedure of quantization was 
correct. Final test would be an agreement of theoretical predictions with real world, i.e. 
experimental/observational data (when the data become available). An agreement at this 
stage would end the story. 

Standard LQC and LQG means the scheme D (DS), nonstandard LQC (our method) 
corresponds to the scheme R (RS). All results (with small exceptions) have been obtained 
so far within DS: almost 25-years of efforts (hundreds of papers) in the case of LQG, and 
almost 10-years of efforts (hundreds of papers) in the case of LQC. There is almost no 
criticism of these results (with small exception: hep-th /0601129, hep-th /0501114 [). 

The above strategy underlies our paper: we intend to verify/confirm the DS LQC results. 
On the other hand, an alternative method may broaden our understanding of various con- 
ceptual issues like, e.g., quantum evolution of a system with the Hamiltonian constraint or 
defining physical observables. Now, let us be more specific. The technical procedure (with 
different interpretation) leading up to Eq. (presenting a Hamiltonian parametrized by 
A) is identical in both DS and RS. The real difference arises when one begins with the 
implementation of the Hamiltonian constraint. 

In DS, one promotes Eq. (|5]) to an operator equation at the quantum level (see, 
gr-qc/0304074, gr-qc/0604013; ff-qc/0607039 ) with the parameter A different from zero^. 
The question arises: Why is the parameter A kept different from zero? The technical answer 
is: In the limit A — )■ the ref^w/ar constraint equation turns into the smf/u/arWheeler-DeWitt 
equation. Users of DS propose the following justification for keeping A 7^ 0: In the loop rep- 
resentation used to quantize the kinematical level (first step of DS) this representation does 
not exist for A = 0. In the second step of DS, one turns the classical constraint equation 
into an operator equation, which one defines consequently in the loop representation of the 
kinematical level. This is why one keeps A 7^ in this operator equation. This procedure 
provokes the next question: Why one uses such an exotic representation which is not defined 
for A 7^ 0? The explanation of the users of DS is the following: such representation is an 



For simplicity, we ignore the difference between tlie so called 'old' and 'improved' schemes notations and 
denote the regularization parameter by the same symbol A. 
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analog of the representation used at the kinematical level of loop quantum gravity (LQG). Is 
it a satisfactory justification? Well, LQG has not been constructed yet: the representation 
of the constraints algebra, based on the achievements of kinematical level, has not been 
found. The problem is extremely difficult because the algebra is not a Lie, but a Poisson 
algebra (one structure function is not a constant but a function on phase space). The users 
of the DS method deeply believe that sooner or later this problem of LQG will be solved and 
LQG will be derived from LQG, so making use of an analogy to LQG is a healthy approach. 
This is a highly promising development, but far from being completed. 

In the RS case, the point of departure (as we said above) is Eq. ([5]). The Hamiltonian 
must be the same, in DS and RS, otherwise the comparison of both methods would be 
impossible. Let us emphasize that Eq. ([5]) is not an effective semi-classical Hamiltonian in 
the RS method. It is considered to be a purely classical, but modified Hamiltonian. One 
can treat it as a one-parameter family of classical Hamiltonians, including the usual general 
relativity Hamiltonian ClS db special case for A = 0. 

It can be seen easily that the singularity becomes 'resolved' at the classical level, for 
A 7^ 0, due to the functional form of Eq. ([S]): big bang turns into big bounce. Why should 
we quantize a cosmological model which is free from the cosmological singularity? 
There are at least three reasons: 

(i) to make comparison with DS results, we must have a quantum model; 

(ii) the parameter A specifying the modification is a free parameter in RS. As the result, 
the classical critical density of matter at the bounce becomes unspecified as it depends on 
A (see, Eq. (ITS]) ). Since it may become arbitrarily big for small enough A, the system may 
enter an arbitrarily small length scale, where quantum effects cannot be ignored. We have 
discussed this issue with all details in [l^ (subsection A of section V); 

(iii) making predictions of our model for quantum cosmic data may be used to fix the free 
parameter A, after such data become available. 

Since the quantum energy density operator turns out to have eigenvalues coinciding with 
the classical expression for the energy density (see, Eqs. f|T3|) and fl22l) ). no new opportunity 
exists to fix A in the context of energy density. Fortunately, the quantum volume operator 
has discrete spectrum depending explicitly on A. Since this results implies that the space 
may have a foamy like structure, the determination of A may be possible (see. Sec. IV). 
Turning the argument around, the detection of such effects could be used to argue that 
quantum gravity effects do occur near the cosmological singularity. 

The RS method reveals that the choice of the exotic Bohr representation at the quantum 
kinematical level of DS in effect corresponds to the setting A 7^ at the classical level. 
After making standard quantization we keep A 7^ because as A — )■ the maximum energy 
density blows up due to Eq. (123|) ). and we get a singularity. Also the spectrum of the 
volume operator becomes continuous, since A := SirG'yhX goes to zero as A — )• 0. Thus, 
any volume of space may smoothly collapse to a singularity. In both RS and DS methods, 
A cannot be treated as a regularization parameter that one should remove from final results. 
Quite the opposite, its presence is necessary for the resolution of the singularity problem. 

Our results agree with the results obtained within the DS method in the sense that the 
classical initial singularity is replaced with the big bounce transition. 

In fact, we have done much more. We have shown that the spectrum of the physical 
volume operator is discrete. Physical means that the spectrum may be used for comparison 
with observational data, in contrast to the case of the spectrum of the kinematical operators. 
We have also made some preliminary steps towards the solution of the evolution problem. 
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which is fundamental for systems with Hamiltonian being the constraint of the system. 

We have found that LQC should be directly linked with experimental data to become 
specific. It is so because A, parametrizing the loop geometry, is a free parameter of LQC. It 
seems, there is no satisfactory way to determine it theoretically within LQC, but forthcoming 
cosmic observations may solve the problem. We have already addressed this issue in the 
context of the DS (see, Sec 4 of Id. Present paper confirms this result within the RS 
approach. Some preliminary agreement with our results (without reference to ours) may be 
found in an updated version of the paper on the robustness of LQC: see, arXiv: 0710.3565, 
version 4, Sec VI. One discusses there 'parachuting by hand' results from full LQC into 
LQC, since the derivation of LQC from LQG has not been obtained yet. If we suitably fix 
A in our method 'by hand', we will get the numerical results concerning the big bounce, 
like the critical energy density of matter field, identical to the ones obtained within the DS 
method (see, Subsec A of Sec V in [12]). 

Our method relies on a direct link with observational data due to the unknown value 
of A. This may make our model suitable for describing observational data despite the fact 
that FRW may have too much symmetry to be a realistic model of the Universe. Taking 
theoretically determined A from incomplete LQG, in DS method, seems to give a model of 
the Universe less realistic than ours. Lacking of theoretically determined numerical value of 
A is not drawback, but advantage of the RS method. 
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